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GAMM  CONFERENCE  ON  NUMERICAL  METHODS  IN 
FLUID  MECHANICS 


Introduction 

The  sixth  Gesellschaft  fiir  ange- 
wandte  Mathematik  und  Mechanik  (GAMM) 
Conference  on  Numerical  Methods  in  Fluid 
Mechanics  was  held  from  25  through 
27  September  1985  in  Gottingen,  West 
Germany.  In  attendance  were  approximate¬ 
ly  200  scientists  from  Europe  as  well  as 
Russia,  India,  Japan,  and  the  US.  The 
meeting  was  hosted  by  Deutsche  For- 
schungs-  und  Versuchsanstalt  fiir  Luft- 
und  Raumfahrt  (DFVLR) .  Papers  dealing 
with  methods  (in  contrast  with  applica¬ 
tions)  were  stressed  in  this  meeting. 
During  the  meeting,  papers  on  Galerkin, 
finite  difference,  finite  volume,  and 
hybrid  methods  for  the  solution  of  the 
Navier-Stokes  and  Euler  equations  were 
presented.  Related  to  these  matters, 
several  papers  dealing  with  grid  genera¬ 
tion,  numerical  dissipation,  and  the 
implementation  of  downstream  boundary 
conditions  were  also  presented.  I  will 
not  be  able  to  report  on  all  of  the 
papers  which  were  presented  because 
there  were  two  parallel  sessions  in  pro¬ 
gress  during  approximately  one-half  of 
the  meeting.  A  proceedings  volume  con¬ 
sisting  of  all  papers  presented  at  the 
conference  will  be  published  by  Vieweg- 
Verlag  in  Xotcs  on  Nxcnevical  Fluid 
.''hanioe  and  should  be  available 
short ly . 


Aero-Acoustic  Modeling 

The  first  paper  was  by  J.C.  Fran- 
jaud  and  J.M.  Hervouet  from  Electricite 
de  France.  Their  interest  was  in  si¬ 


lencing  piping  systems  for  carrying 
steam  from  the  boilers  to  the  generators 
in  steam  power  plants.  They  needed  a 
way  to  predict  the  aerodynamic  noise 
produced  in  such  systems.  Initially, 
they  set  up  a  test  in  a  transparent  test 
rig  in  which  hot-wire,  acoustic,  laser 
Doppler  velocimetry  (LDV) ,  and  inter¬ 
ferometric  measurements  were  made.  The 
working  fluid  was  air.  Various  compo¬ 
nents  present  in  steam  distribution  sys¬ 
tems  could  be  inserted  upstream  of  the 
test  section.  In  the  case  of  an  orifice. 


they  found  that  the  flow  leaving  the 
orifice  reattached  to  the  upper  side  of 
the  test  section,  producing  a  highly 
three-dimensional  asymmetric  condition. 

In  an  attempt  to  simulate  this 
test,  they  used  a  hydrodynamic  code 
which  they  modified  for  compressibility 
effects.  To  simulate  the  turbulence, 
the  k-e  turbulence  model  was  used.  A 
fractionary  step  integration  approach 
was  used  in  which  the  advectlon  terras 
were  Integrated  by  a  characteristic 
method.  Initially,  their  calculations 
were  carried  out  using  a  law-of-the-wall 
model.  However,  this  model  gave  a  sym¬ 
metric  rather  than  the  desired  asymme¬ 
tric  solution.  When  they  refined  the 
grid,  discarded  the  law-of-the  wall 
model,  and  carried  out  their  calcula¬ 
tions  all  the  way  to  the  wall,  they 
found  the  asymmetric  results  which  they 
had  seen  in  their  experiments.  For  both 
calculations  they  used  the  measured  up¬ 
stream  velocity  profile  as  inlet  bound¬ 
ary  conditions.  This  Illustrates  the 
Importance  of  correctly  modeling  the 
wall  boundary  conditions. 

A  second  calculation  was  made  which 
Involved  the  flow  through  a  dual-throat 
nozzle.  A  two-dimensional  Euler  code 
was  used  which  employed  a  fractionary- 
step  integration  approach,  similar  to 
that  used  in  the  hydrodynamics  code.  In 
order  to  accurately  predict  the  acoustic 
features  of  the  flow,  a  new  turbulence 
model  was  developed.  It  Included,  in 
addition  to  the  standard  turbulence  ki¬ 
netic  energy  and  turbulence  dissipation 
quantities,  transport  equations  for  two 
new  variables:  the  acoustic  correlation 
and  an  energetic  variable.  Although  a 
new  acoust  Ic-oriented  turt)ulenrp  n:odel 
was  developed,  no  results  were  availnblo 
to  validate  its  performance.  Neverthe¬ 
less,  some  careful  thought  seems  to  have 
gone  into  the  constnu-tion  of  their  tur¬ 
bulence  model.  Because  of  this,  further 
developments  of  this  research  should  he 
closely  watched. 

Transit  ion  Control 

The  second  paper  was  given  hv  F. 
l.aurien,  of  the  DFVl.R,  Gottingen.  He 
described  his  spectral  calculations  of 


the  transition  mechanisms  responsible 
for  the  laninarizing  effects  of  suction 
and  blowing.  In  order  to  produce  the 
onset  of  turbulence,  he  introduced  spa¬ 
tially  periodic  disturbances.  The  evolu¬ 
tion  of  these  disturbances  was  modeled 
with  the  three-dimensional,  time-depend¬ 
ent,  I.'avier-Stokes  equations  in  a  coor- 
d.inate  system  which  moved  downstream  at 
the  speed  of  the  Tollmien-Schlichting 
waves.  U’hat  his  calculations  showed  was 
that  blowing  and  suction  produced  a  sig¬ 
nificant  decrease  in  the  turbulence 
fluctuations  compared  with  Che  growth 
predicted  when  no  boundary-layer  control 
was  present.  In  addition,  the  simulated 
effects  of  blowing  and  suction  on  the 
production  of  turbulent  fluctuations 
agreed  well  with  experimental  measure¬ 
ments.  Both  the  well-known  harmonic  and 
subharmonic  transition  structures  were 
predicted  by  Che  code  and  it  was  possi¬ 
ble  to  extract  from  Che  model  the  de¬ 
tails  of  the  flow  mechanisms  responsi¬ 
ble  for  Che  transition  phenomenon.  The 
calculations  showed  chat  transition  can 
’e  delayed  if  suction  or  blowing  is  pro¬ 
vided  early  enough  in  Che  transition 
process.  This  was  substantiated  by  shear 
naps,  which  were  clearly  devoid  of  the 
tongi'.e-Mke,  three-dimensional  struc¬ 
tures  which  are  known  precursors  of  tur¬ 
bulence.  This  is  not  to  say  that  the 
ctfert  of  the  disturbances  had  been  re¬ 
moved,  Span-wise  vortices  were  still 
present;  however,  they  did  not  develop 
into  the  tvpes  of  structures  which  are 
knowai  to  lead  to  turbulence.  This  indi¬ 
cates  that  transition  control,  while 
deliving  transition,  certainly  cannot 
return  the  *■  ’  ow  to  Its  initial  uniform 

s  t  c.  r  p  . 

Vo •'tex  Shedding 

P.  Kiehm  c'f  the  Institute  for  Ther- 
n’a  I  and  Fluid  Dynamics  at  the  Ruhr  Uni- 
versitv  In  Bochum,  West  Germany,  de- 
scril'e.l  his  use  of  3  marker  and  cell 
(MAGI  technique  to  calculate  the  flow 
around  a  confined  circular  cylinder. 
G<  n  ‘  Igurat  i '.’OS  where  the  cylinder  was 
bor',  .a’nllned  by  planes  perpendicular 
and  oaruillel  to  the  axis  of  the  cylinder 
were  exam’red.  ’'pstrear  of  the  cylinder 


the  velocity  profile  was  assumed  to  be 
parabolic.  Of  particular  interest  was 
the  Reynolds  number  at  which  periodic 
vortex  shedding  was  observed.  It  was 
found  that  confining  the  cylinder  pro¬ 
duced  a  significant  increase  in  the  Rey¬ 
nolds  number  at  which  periodic  shedding 
was  observed;  for  high  blockage  ratios, 
Reynolds  numbers  as  high  as  67  were 
found  compared  with  Reynolds  numbers  of 
40  for  an  unconfined  cylinder.  Also, 
confinement  was  found  to  damp  out  the 
wake  flow  downstream  of  the  cylinder  and 
return  the  velocity  profiles  to  their 
initial  parabolic  form  more  quickly  than 
in  the  unconfined  case.  In  addition, 
for  the  axially  confined  cylinder,  large 
velocity  components  In  the  direction  of 
the  cylinder  axis  were  predicted.  This 
contributed  to  a  strong  helical  motion 
which  occurred  downstream  of  the  cylin¬ 
der,  see  Figure  J . 

Domain  Decomposition  Methods 

W.  Fritz  (Dornier)  talked  about  his 
multiblock  approach  for  calculating  two- 
dimensional  flow  around  complex  geome¬ 
tries.  The  concept  here  was  to  subdivide 
the  flow  field  into  a  number  of  topolog¬ 
ically  simple  subregions  to  simplify  the 
grid  generation  process.  Both  invlscid 
and  boundary-layer-corrected  invlscid 
calculations  were  shown.  One  of  the 
calculations  which  he  showed  was  the 


Figure  1.  Confined  circular  cylinder. 
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flow  over  an  automobile.  In  this  case 
the  Reynolds  number  ranged  from  1x10^  to 
5x10®.  Both  the  viscous  and  boundary- 
layer-corrected  calculations  predicted 
separated  flow  at  the  rear  of  the  body. 
The  phenomenon  responsible  for  the  sepa¬ 
ration  was  completely  different  in  the 
two  cases,  however.  In  the  case  of  the 
viscous  calculations,  the  separation  was 
due  to  the  existence  of  an  unfavorable 
pressure  gradient  at  the  rear  of  the 
automobile.  This  resulted  in  the  sepa¬ 
ration  point  being  located  upstream  of 
the  base  of  the  body  and  produced  a 
value  of  the  base  pressure  which  was  in 
reasonable  agreement  with  experiment. 
In  the  invlscld  calculations,  on  the 
other  hand,  the  zone  of  separation  was 
restricted  to  the  base  region  Itself  and 
was  due  to  the  sharp  corners  at  the  edge 
of  the  base  region.  These  corners  pro¬ 
duced  a  total  pressure  loss  (even  though 
the  calculations  were  invlscid)  because 
of  the  discretization  errors  produced  by 
the  rapidly  changing  metrics.  Not 
unexpectedly,  the  base  pressure  predica¬ 
ted  by  the  inviscld  calculations  was  in 
poor  agreement  with  experiment. 

By  far  the  most  sophisticated  ap¬ 
plication  of  the  multlblock  approach  was 
presented  by  E.  Venkatapathy  (NASA, 
Ames).  Instead  of  the  term  multiblock 
he  preferred  to  use  "patched  grids,"  and 
hr  defined  several  ways  in  which  patch¬ 
ing  could  take  place.  The  advantages  of 
patched  (multiblock)  grids  are  that  they 
allow  good  alignment  between  the  grid 
and  the  flow  and  a  fine  mesh  can  be  pro¬ 
duced  exactly  where  it  is  needed.  On 
the  other  hand,  difficulties  are  intro¬ 
duced  by  this  approach,  primarily  be¬ 
cause  of  the  increased  complexity  of  the 
data  structure.  This  is  because  with 
overlapped  grids  there  is  a  need  for 
both  global  and  local  indexing  and  the 
computed  data  must  be  passed  between  the 
global  and  local  meshes  in  an  efficient 
fashion.  The  numerical  algorithm  which 
was  used  In  the  calculations  was  the 
conservative  supra  characteristic  method 
(CSCM) ,  which  is  a  diagonally  dominant, 
approximate-factorization  technique.  It 
has  the  Important  characteristic  that  it 
is  an  upwind  scheme.  For  patching  meth¬ 


ods,  upwinding  is  an  absolute  require¬ 
ment.  Linear  interpolation  was  used  to 
exchange  data  between  the  grid  systems. 

The  fine  mesh  of  the  local  grid  has 
less  numerical  diffusion  and  thus  pro¬ 
duces  much  sharper  shock  calculations. 
Venkatapathy  showed  a  calculation  of 
Mach  5.0  flow  over  a  cone-cylinder  which 
clearly  demonstrated  the  improvements 
which  can  be  obtained  with  patched 
grids.  A  comparison  was  shown  between  a 
coarse  (26x26)  grid  on  which  patching 
was  employed  and  a  fine  grid  (101x101) 
in  which  no  patching  was  used.  The  re¬ 
sults  with  patching  on  the  coarse  grid 
were  at  least  as  good  as  the  fine  mesh 
calculations.  The  patching  technique 
can  be  automated,  thus  producing  an 
evolutionary,  adaptive  grid  structure 
which  adjusts  to  the  complexities  of  the 
flow  as  they  are  encountered  in  the  cal¬ 
culations. 

A  particularly  interesting  applica¬ 
tion  of  the  patching  approach  was  pre¬ 
sented  by  M.  Ral  of  the  Information  Gen¬ 
eral  Corporation,  Mountain  View,  Cali¬ 
fornia.  He  showed  a  rotor-stator  inter¬ 
action  problem  in  which  the  rotor  grid 
moved  relative  to  the  stator  grid.  In 
his  presentation,  he  highlighted  the 
treatment  of  the  interface  points.  He 
Indicated  that  interface  procedures  must 
be  stable,  accurate  in  space  and  time, 
conservative,  applicable  to  generalized 
coordinate  systems,  and  easy  to  incor¬ 
porate  in  existing  codes.  For  geometric 
simplicity,  he  used  circular-arc,  bicon¬ 
vex  airfoils  for  the  rotor  and  stator 
blades.  Identical  H-grlds  were  used  for 
both  the  stationary  (stator)  and  the 
moving  (rotor)  patches.  Osher's  total 
variation  diminishing  (TVD)  scheme  was 
used  to  solve  the  unsteady  Euler  equa¬ 
tions.  A  spectacular  color  motion  pic¬ 
ture  showing  the  case  of  a  1.5  Mach  num¬ 
ber  inlet  flow  Illustrated  the  results 
of  his  calculations.  The  interaction 
between  the  rotor  and  the  stator  flow 
fields  was  dramatically  evident.  The 
resolution  of  the  calculations  seemed 
excellent,  particularly  in  the  vicinity 
of  the  shock-shock  Interactions.  A  seg¬ 
ment  showing  a  "zoomed-ln"  view  of  the 
flow  at  the  leading  edge  of  the  stator 


plates  clearly  showed  the  periodic  at¬ 
tachment  and  the  detachment  of  the  shock 
in  response  to  the  flow  generated  by  the 
rotor.  There  was  no  visible  evidence  of 
the  patching  interface  in  any  of  the 
calculations.  In  addition  to  turboma¬ 
chinery  problems,  this  method  can  be 
applied  to  helicopter  blade  flows,  prop- 
fan  calculations,  and  propeller/fuselage 
interaction  problems. 

TVD  Calculations 

H.C.  Yee  (NASA,  Ames)  presented  a 
new  symmetric-differencing,  total  varia¬ 
tion  diminishing  (TVD)  scheme  and  com¬ 

pared  it  with  the  original  upwind  TVD 
scheme  developed  by  Osher.  Her  calcu¬ 
lations  showed  that  a  30  to  40  percent 
savings  in  computational  time  was  possi¬ 
ble  with  the  new  method. 

In  another  TVD  calculation,  C.M. 
Kwong  (University  of  Salford,  UK)  high¬ 
lighted  an  Important  advantage  of  TVD 
methods;  namely,  that  shocks  can  be 
captured  without  the  spurious  oscilla¬ 

tions  which  often  occur  when  the  Mac- 
Cormack  algorithm  is  used.  Kwong  showed 
calculations  of  flow  in  a  shock  tube, 
flow  in  a  channel  with  a  bump,  and  an 
oblique  shock  reflection  calculation  in 
which  such  oscillations  were  notably 
absent.  Because  TVD  schemes  offer  a 

rational  method  for  controlling  the 

proper  amount  of  dissipation,  no  prob¬ 
lem-dependent  adjustments  were  needed  to 
achieve  this  result. 

Viscous/Inviscid  Interaction 

M.  Schmatz,  Messerschmitt-Bolkow- 
I'lohm  (MBB)  ,  presented  his  strong  vis¬ 
cous /ir.viscld  interaction  airfoil  cal¬ 
culations.  His  approach  was  to  divide 
the  flow  field  into  three  zones;  a  zone 
upstream  but  far  away  from  the  airfoil 
where  the  fuler  equatiiMis  are  used,  a 
regirn  upstream  but  near  the  airfoil 
where  a  disp! acement-thickness-corrected 
'^iilcr  calculation  is  used,  and  a  region 
ihjwnstream  of  the  airfoil  where  the 
Navler-Stokes  equatic'ns  are  used.  This 
is  shc'wn  in  Figure  2.  The  equivalent 
Invisc'd  ‘'low  which  results  from  the 
’'uler/bound.irv-layer  calculations  is 
iter.itivelv  coupled  with  the  Navler- 


Flgure  2.  Overlapping  subdomains. 


Stokes  solution  by  means  of  a  modified 
Schmatz  alternating  algorithm.  A  Bald- 
wln-Lomax  algebraic  two-layer  stress 
model  was  used  for  the  turbulence  cal¬ 
culations.  Results  were  shown  at  Mach 
numbers  of  0.504,  0.72,  and  0.83  at 
three  angles  of  attack  and  at  Reynolds 
numbers  of  approximately  6x10®.  Some 
difficulties  were  encountered  in  accur¬ 
ately  predicting  the  skin  friction,  and 
wiggles  appeared  downstream  of  the 
region  where  the  Euler  and  Navler-Stokes 
codes  were  coupled.  To  remove  these 
wiggles  it  was  suggested  that  an  impli¬ 
cit  Euler  code  or  characteristic-orien¬ 
ted  Navier-Stokes  code  could  be  used. 
The  big  advantage  of  Schmatz 's  method 
lies  in  three-dimensional  calculations 
and  not  in  the  two-dimensional  cases 
which  were  presented  here.  Future  plans 
include  extending  this  work  to  two- 
dimensional  separated  flows  and  then  to 
three-dimensional  calculations. 

D.  Chaussee  (NASA,  Ames)  presented 
his  parabolized  Navier-Stokes  calcula¬ 
tion  for  the  flow  over  a  cone.  A  Bald- 
win-Lomax  turbulence  model  was  used. 
Shock  waves  were  fitted  using  charac¬ 
teristic  theory  to  reduce  the  number  of 
points  needed  in  the  calculations.  Cal¬ 
culations  at  4-,  10-,  and  14-degree 
angles  of  attack  were  made  with  as  many 
as  121x45  points  in  the  cross  plane. 
Even  with  the  finest  grid  it  was  impos¬ 
sible  to  obtain  grid-insensitive  solu¬ 
tions,  especially  on  the  pressure  sur¬ 
face  of  the  cone.  Accordingly,  a  run 
with  in  points  was  made.  This  calcula¬ 
tion  required  approximately  2.2  hours  on 
a  CRAY  X-MP. 


Euler  Solutions 

L.  Eriksson  of  the  Aeronautical  Re¬ 
search  Institute  of  Sweden  (FFA)  showed 
how  it  was  possible  with  local  grid  re¬ 
finement  to  greatly  increase  the  accur¬ 
acy  of  Che  calculation  without  greatly 
increasing  the  number  of  computational 
points.  First,  a  coarse  mesh  solution 
was  obtained,  then,  in  the  region  selec¬ 
ted  for  local  grid  refinement,  values  of 
the  flow  variables  were  estimated  either 
by  linear  means  or  by  means  of  extra¬ 
polation  of  the  appropriate  Rlemann  in¬ 
variants.  In  a  calculation  made  on  a 
N'ASA  0012  airfoil,  25  percent  more 
points  were  added  in  the  vicinity  of  the 
nose.  The  surprising  result  was  that 
this  reduced  Che  total  pressure  error  by 
75  percent.  In  addition,  grid  refine¬ 
ment  produced  a  sharpened  shock  struc¬ 
ture  which  in  some  cases  persisted  even 
beyond  the  zone  selected  for  grid  re- 
f inement . 

A.  Rizzi  CFFA)  described  his  invis- 
cid  calculations  of  the  flow  over  a 
cranked,  twisted,  delta  wing.  The  pur¬ 
pose  of  his  calculations  was  to  study 
unstable  flows.  Such  calculations  are 
onlv  possible  if  there  are  modes  in  such 
flows  which  do  not  experience  explosive 
growth.  For  these  calculations  an  expli¬ 
cit,  finite-volume  technique  was  used 
employing  second-  and  fourth-order  arti¬ 
ficial  dissipation.  Runge-Kutta  integra¬ 
tion  was  used.  Calculations  were  pre¬ 
sented  for  a  Mach  number  of  1.2  at  an 
angle  of  attack  of  10  degrees.  Both 
medium  mesh  (80,000  points)  and  fine 
mesh  (540,000  points)  calculations  were 
made;  however,  it  was  not  possible  to 
obtain  a  mesh-independent  solution, 
based  on  the  level  of  artificial  viscos¬ 
ity  which  was  present  in  the  flow,  the 
probable  value  of  the  Reynolds  number 
was  approximately  1,000. 

Summary 

At  the  end  of  the  meeting  the  ses¬ 
sion  chairmen  were  invited  to  give  their 
reactions  to  the  papers  which  had  been 
presented  and  to  summarize  for  the  audi¬ 
ence  their  perceptions  of  the  state  of 


the  art  in  computational  fluid  dynamics. 
They  generally  agreed  that  three-dimen¬ 
sional  Euler  simulations  involving  up¬ 
wards  of  250,000  points  were  common¬ 
place,  but  whether  or  not  these  codes 
presently  exist  in  versions  which  could 
be  used  for  production  calculations  is  a 
matter  of  opinion.  However,  within  the 
next  decade  the  necessary  degree  of  user 
friendliness  for  production  applications 
will  certainly  be  developed.  The  outlook 
for  Navier-Stokes  codes  is  less  clear. 
Production  codes  employing  simple  turbu¬ 
lence  models  should  be  available  within 
the  next  decade.  These  codes  will  enable 
three-dimensional  calculations  involving 
more  than  eight  million  points  to  be 
accomplished  and  will  most  certainly 
make  use  of  large  eddy  simulations. 

Other  chairmen  expressed  a  more 
conservative  viewpoint  and  stressed  the 
need  for  new  algorithms  (for  both  vector 
and  parallel  computations)  for  inversion 
of  large  matrices.  They  expressed  the 
need  for  more  accurate  Euler  schemes  to 
minimize  the  effects  of  numerical  dissi¬ 
pation.  They  felt  that,  for  some  time 
to  come,  multizone  decomposition  employ¬ 
ing  some  type  of  adaptive  scheme  will  be 
an  alternative  to  Navier-Stokes  calcula¬ 
tions  for  practical  computations. 

Whatever  the  code,  attention  must 
be  given  to  reducing  its  development 
time.  There  must  be  better  reporting  of 
one  another's  mistakes  and  of  the 
strength  and  weaknesses  of  methods  now 
in  existence. 

A  final  comment  was  made  bv  M. 
Napolitano  of  the  University  of  Rarl , 
Italy,  who  said  that  although  inflation¬ 
ary  claims  of  code  capability  needed  to 
be  avoided,  inflationary  code  utiliza¬ 
tion  should  not.  "You  need  to  be  criti¬ 
cal  of  a  code"  he  said,  "but  not  too 
cautious."  Often,  he  said,  a  new  ap¬ 
proach  Is  discovered  by  examining  the 
failures  of  a  calculation  which  is 
pushed  too  far.  He  cited,  as  an  exam¬ 
ple,  the  failure  of  conventional  bound¬ 
ary  layer  techniques  which  when  applied 
to  separated  flows  led  to  the  develop¬ 
ment  of  indirect  boundarv-1 aver  methods. 


